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If all the supersymmetry particles (sparticles) except a light Higgs boson are too heavy to be 
directly produced at the Large Hadron Collider (LHC) and Tevatron, a possible way to reveal 
evidence for supersymmetry is through their virtual effects in other processes. We examine such 
supersymmetric QCD effects in bottom pair production associated with a light Higgs boson at the 
LHC and Tevatron. We find that if the relevant sparticles (gluinos and squarks) are well above 
the TeV scale, too heavy to be directly produced, they can still have sizable virtual effects in this 
process. For large tan/3, such residual effects can alter the production rate by as much as 40 percent, 
which should be observable in future measurements of this process. 
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D I. INTRODUCTION 

Searching for the Higgs boson is one of the most important tasks in particle physics. The existence of a relatively 
light Higgs is suggested by high precision fits to the data in the Standard Model (SM) and also is theoretically favored 
in the Minimal Supersymmetric Standard Model (MSSM) [1]. Verification of the existence of a light Higgs boson 
at the LHC or Tevatron is therefore a very important test for both the SM and the MSSM. Among the various 
production channels for a light Higgs boson at the hadron colliders, the production in association with a bottom 
quark pair, pp(pv pp) — > hbb + X, plays an important role in testing the bottom quark Yukawa couplings. While 
this process has a small cross section in the SM, in the MSSM this production mechanism can be a significant source 
of Higgs bosons since the bottom quark Yukawa coupling in the MSSM is proportional to tan/3 (defined as i^/fi 
with ui,2 being the vacuum expectation values of the two Higgs doublets) and the current analyses favor large tan/3. 
This process has been studied at next-to-leading order in perturbative QCD [2-8]. Due to the importance of this 
production mechanism in testing the bottom quark Yukawa coupling in the MSSM, the supersymmetry (SUSY) loop 
effects in this process should also be considered. Of course, among the SUSY loop effects the one-loop SUSY QCD 
corrections are the dominant. 

In the present work we examine the one-loop SUSY QCD effects in this process. Instead of performing a complete 
one-loop calculation, which would be quite complicated since it involves many five-point box diagrams, we focus on 
the so-called SUSY residual effects, i.e., the SUSY effects in the heavy limit (<; TeV) of the sparticles involved. Our 
motivations are the following: 

• It is possible that all the sparticles except a light Higgs boson are too heavy to be directly produced at the 
LHC and Tevatron, such as in the split SUSY scenario proposed recently by Arkani-Hamed and Dimopoulos [9]. 
Although the fermionic partners (gauginos and Higgsinos) in this scenario are required to be relatively light in 
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order to ensure gauge coupling unification and provide dark matter candidates, they are not necessarily below 
a TeV, as recently shown [10]. Thus, it is possible for the LHC and Tevatron to observe no sparticles except a 
light Higgs boson. In that case a possible way to reveal a hint of supersymmetry is through its residual effects 
in observable processes. 

• Unfortunately (or fortunately), SUSY virtual effects decouple in most processes when SUSY particles become 
very heavy. However, we know that the only processes where SUSY has residual effects are those processes 
involving Higgs-fermion Yukawa couplings, as first studied for h — > bb decay [11,12], and also for certain pro- 
duction processes [13,14]. The production reaction pp(ov pp) — > hbb + X at the LHC or Tevatron is well suited 
for revealing SUSY residual effects since it involves the hbb coupling. Compared with pp(or pp) — > hb + X, this 
process is also easier to detect since it contains an extra hard b jet in the final state. Once the Higgs boson (h) 
is observed and its mass is measured through other processes such as gluon-gluon fusion, this reaction can be 
used to measure the bottom quark Yukawa coupling and to observe the expected residual effects of SUSY. 

Note that the existence of SUSY residual effects in some Higgs process does not mean SUSY is not decoupling in low 
energy processes. As shown in previous studies, and also in the following work, the residual effects exist when Ma 
remains light; when Ma is heavy, together with all other SUSY masses, the residual effects do vanish. 

This paper is organized as follows: In Section II we present our strategy for the calculation of the one-loop SUSY 
QCD corrections. In Section III we perform numerical calculations and obtain the residual effects in the limit of heavy 
SUSY masses. The conclusion is given in Section IV and the detailed analytic formulas obtained in our calculations 
are presented in the Appendix. 



II. ONE-LOOP SUSY QCD CORRECTIONS 

The production pp(ov pp) — > hbb + X proceeds through the parton-level processes gg — > hbb and qq — > hbb. The 
one-loop SUSY QCD corrections to this process have a huge number of one-loop diagrams, including many five-point 
box diagrams. However, among all these diagrams only the one-loop diagrams involving the bottom quark Yukawa 
coupling have residual effects as the SUSY masses become very heavy. Therefore, in our calculation we need only 
consider the loop corrected bottom quark Yukawa coupling diagrams shown in Fig. 1 and Fig. 2. 

In our loop calculations we used dimensional regularization to control the ultraviolet divergences and adopted the 
on-mass-shell renormalization scheme. Each effective hbb vertex in Fig.l contains two parts: one is the irreducible 
three-point vertex loop contributions and the other is the counterterms 8V hb i — g^^SZ with g^^ denoting the tree-level 
hbb coupling and SZ is the renormalization constant given by 



SZ L SZ R 5m b 



(1) 



where SZl^r and Snib are respectively the renormalization constant for the b quark field and mass. They can be 
extracted from the one- loop self-energies shown in Fig. 2(b) by using the on-mass-shell renormalization condition. 
They are given by 
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where Aj = af + b\ , Af = of — b\ and A\ = 2dibi with aj and bi given in the Appendix. B ,i are the 2-point Feynman 
integrals given in [15], and their functional dependence is indicated in the brackets following them. 

The counterterm is universal for each hbb vertex shown in Fig. 1. However, although the irreducible three-point 
vertex loops have the same topological structure shown in Fig. 2(a), the results are different for different hbb vertices 
in Fig. 1 because they depend on the external momenta. The results are lengthy and are presented in the Appendix. 
We have checked that all the ultraviolet divergences do cancel as a result of renormalizability of the MSSM. 

Note that in our calculations we adopted the so-called on-mass-shell scheme, in which the renormalized mass 
ra\, = m\ — Snib (mjj is the bare mass) is the physical mass, i.e., the pole of the 6-quark propagator [11,16,17]. The 
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FIG. 1. Feynman diagrams for parton-level subprocesses for pp(px pp) — > hbb + X with one-loop SUSY QCD corrected hbb 
vertices. The large dots denote the effective hbb vertex with one-loop SUSY QCD corrections. 



main difference of this scheme with the MS scheme is that in the MS scheme a running b-quark mass is introduced to 
absorb the leading part of the corrections (for example, the large logarithms in QCD corrections) [16]. In calculating 
SUSY-QCD corrections, there are no such large logarithms and the on-mass-shell scheme is usually adopted ( an 
extensive discussion about this issue was provided in [17]) . 

Including the one-loop SUSY-QCD corrections to the bottom quark Yukawa coupling, the renormalized amplitude 
for pp(ov pp) — > hbb + X can be written as 

M = Mf + SM™ + Ml* 9 + 8M" , (5) 

where Mq and 5M represent the tree-level amplitude and one-loop SUSY-QCD corrections, respectively. The detailed 
expression for SM is given in the Appendix. 

In our calculation we performed Monte Carlo integration to obtain the hadronic cross section by using the CTEQ5L 



1 If we use the MS scheme in calculating SUSY-QCD corrections to the hbb coupling, i.e., define a running b-quark mass to 
absorb some SUSY-QCD correction effects, we will obtain the approximately same result for the ratio of the cross sections with 
and without SUSY-QCD corrections ( because by doing this we merely relocated some correction effects and the total one-loop 
SUSY-QCD correction effects are not changed). 
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(a) (b) 

FIG. 2. Feynman diagrams for the one-loop SUSY QCD corrections to the hbb vertex: (a) the irreducible vertex loops ana 
(b)the self-energy loops of the b guark. 



parton distribution functions [18] with Q = and requiring the transverse momentum of the two b-jets to be larger 
than 15 GeV. 

To exhibit the size of the corrections we define the ratio 

Asqcd = , (6) 

where oo is the tree-level cross section. 

III. NUMERICAL RESULTS 

Before performing numerical calculations, we must choose the parameters involved. For the SM parameters we used 
m w = 80.448 GeV, m z = 91.187 GeV, m t = 178 GeV, m b = 4.5 GeV, sin 2 8 W = 0.223, and the two-loop running 
coupling constant a s (Q). For the SUSY parameters, apart from the charged Higgs mass, gluino mass and tan/3, the 
mass parameters of sbottoms are involved. The mass-squared matrix for the sbottoms takes the form [19] 

„ / to? rribXl \ , 



where 



m L = m % + m l- m2 z < \-\ sin2 °w) cos(2/3) , (8) 



1 

3 

X b = A b - /xtan/3 , (10) 



m l = m D +m b ~ -m 2 z sin 2 9 W cos(2/3) , (9) 



Here mj~ and m 2 ^ are soft-breaking mass terms for the left-handed squark doublet Q and the right-handed down 

squark D, respectively. A b is the coefficient of the trilinear term HiQD in the soft- breaking terms and \i is the 
bilinear coupling of the two Higgs doublets in the superpotential. Thus, the SUSY parameters involved in the 
sbottom mass matrix are TOq, to^, Ab, [i and tan/3. 

The mass-squared matrix is diagonalized by a unitary transformation which relates the weak eigenstates b^ ^ to 
the mass eigenstates &i 2 : 



h \ _ / cos6> fc sin 6* b \ / br, 
b 2 J \ - sin 8b cos Ob ) \b R 



(11) 



with the mixing angle and masses determined by 
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In addition, we consider the following experimental constraints: 

(1) \x > and large tan/3, in the range 5 < tan/3 < 50, which are favored by the recent muon g — 2 measurement [20]. 

(2) The LEP and CDF lower mass bounds on gluino and sbottom [21] 

m bi > 75.0 GeV, m~ g > 190 GeV . (14) 

Next, we present numerical results for LHC (y^i = 14 TeV) and Tevatron (^/s — 2 TeV) in two representative cases: 

Case A: All SUSY mass parameters are of the same size and much heavier than the weak scale, but Ma is 
fixed at weak scale, i.e., 
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FIG. 3. SUSY QCD effects Asqcd = (cr - a )/a versus the common SUSY mass Msusy (= Mq = M^ = A b = M- g = fi) 
for a fixed Ma and different values o/tan/3. The corresponding range of mh is 94 ~ 113 GeV for tan ft = 5, 102 ~ 121 GeV 
for tan/3 = 15 and 103 ~ 122 GeV for tan/3 = 30. 



In this case the mixing of sbottoms is maximal, i.e., 9 b ~ ±7r/4. Figs. 3 and 4 show the dependence on the SUSY 
scale Msusy for different values of tan/3 and Ma- We see that Asqcd approaches a non- vanishing constant as 
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FIG. 4. Same as Fig. 3 but for a fixed value o/tan/3 and different values of Ma- The corresponding range of mu is 99 ~ 118 
GeV for m A = 200 GeV, 100 ~ 119 GeV for m A = 300 GeV and 100 ~ 120 GeV for m A = 500 GeV. 



Msusy becomes large. The effects are enhanced by tan/3. For tan/3 = 30, the residual effects can be as large 
as 40% 2 . Fig. 4 shows that as Ma becomes large the size of residual effects decrease. 

Case B: All SUSY mass parameters, including Ma, are much larger than the weak scale, i.e., 

Msusy = Mq =M f) = A b = M~ g =» = M A > M EW . (16) 

This case also gives maximal mixing for sbottoms. In Fig. 5 we show the dependence on the MSSM mass scale 
Msusy- This figure explicitly exhibits the decoupling behaviour of SUSY QCD as Msusy becomes large. 

The underlying reason for the existence of SUSY residual effects in fermion Yukawa couplings when Ma is fixed 
at the weak scale is that certain couplings are proportional to SUSY mass parameters [22]. As pointed out in [23], 
the SUSY QCD residual effects in the hbb coupling are proportional to tan j3. For hit coupling, however, the SUSY 
QCD residual effects are proportional to cot (3. For the process pp(or pp) — > hit + X the calculation of SUSY QCD 
corrections is analogous to the pp(or pp) — > hbb + X case. However, we found that the SUSY residual effects in 
pp(or pp) hit + X are quite small, reaching only 3% for tan/3 = 5 and being smaller for larger tan/3 values. Such 
a small effect is even less than the next-to-leading-order QCD corrections [24-29] and is not likely to be observed at 
the LHC or Tevatron. 



2 Note that when the one-loop effects are very large, higher order loops should also be considered. In Ref. [17] resummation 
techniques are proposed to improve the one-loop results. 
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FIG. 5. Same as Fig.3 but with Ma = Msusy- The corresponding range of m h is 94 ~ 116 GeV for tan = 5, 102 ~ 122 
GeV for tan/3 = 15 and 103 ~ 123 GeV for tan/3 = 30. 



Note that in our calculations, instead of using the effective hbb vertex [12], we performed the complete one-loop 
calculations for hbb Yukawa coupling, which, as pointed out in the paragraph following eq.(4), are dependent on the 
external momenta of hbb vertex. We found that when SUSY mass scale is larger than about 1 TeV, using the effective 
hbb vertex [12] is a good approximation. 

As we pointed out in the Introduction, we only considered the one-loop SUSY QCD corrections, which are the 
dominant part of the SUSY corrections. Of course, the SUSY electroweak corrections to hbb coupling also have 
residue effects, which, however, are obviously suppressed by a factor aEw/o^s relative to the SUSY QCD effects. 
Since much more parameters are involved in the SUSY electroweak sector, we did not perform a detailed calculation 
for the SUSY electroweak corrections here. 

IV. CONCLUSION 

We examined the supersymmetric QCD effects in bottom pair production associated with a light Higgs boson at the 
LHC and Tevatron. We found that when the relevant sparticles are heavy, well above the TeV scale, they nevertheless 
contribute sizable virtual effects in this process. For large tan/3, these residual effects can alter the production cross 
section by as much as 40 percent and thus should be measurable in observations of this process at the LHC or 
Tevatron. 

If only one light Higgs boson, or perhaps several Higgs bosons, are discovered at the LHC or Tevatron, it could 
indicate that the SUSY scale is quite high; above the TeV scale. Of course, the fine-tuning problem then remains 
and supersymmetry loses one of its merits. But this might possibly happen since supersymmetry does not have to 
solve the fine-tuning problem, as argued recently by Arkani-Hamed and Dimopoulos [9]. In that case, a possible way 
to reveal a hint of supersymmetry is through its residue effects. Higgs production associated with a pair of bottom 
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quarks is then well suited for such a seeking a clue for supersymmetry. 
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APPENDIX A: EXPRESSIONS OF SUSY QCD CORRECTIONS TO THE AMPLITUDE 

The coupling at the vertex h°bib* is needed in our calculations. It is given by 

V(h%b*) = igQij , (Al) 



where 



Qn — Ci cos 2 6b + c 2 sin 2 9b + 2c 3 sin 9b cos 9b , (A2) 

Q12 = (c2 - ci) sin 9 h cos 9 b + c 3 (cos 2 9 b - sin 2 9 b ) , (A3) 

Q21 = (C2 - ci) sin^cos^ + c 3 (cos 2 9 b - sin 2 9 h ) , (A4) 

Q22 = ci sin 2 9b + C2 cos 2 9b — 2c 3 sin 9b cos 9b , (A5) 



with 



ci = — ^f- (±-\sin 2 9 w )sm(a + P) , (A6) 
cos 6^ \2 3 / 

c 2 = ^f- \ sin 2 <9 W sin(a + 0) , (A7) 

COS (/jy 3 

c 3 = (Ab sin a + u cos a) . (A8) 

2mw cos p 

Here a is the mixing angle between the two neutral CP-even Higgs bosons. Then one-loop contributions 8M qq and 
5M" are given by 

SM " q = ~wu^f^ Tkl 2 M Pi)ru( P2 )u(k 2 ) [(a, - 6 j7 5 )(- faCu- foC 12 +m~ g C ) 

S[(K 2 + fc 3 j 771;,] 



x ( a j + b j"f 5 )Qij - G5Z] {Jt 2 + fa + mbhij,v(ki)(-k2,-k3,mg,m~ bi ,rn~ b .) , (A9) 
5M ^ = -77r—T^^v{p 1 ) 1 »u{p 2 )u{k 2 ) 1 ^fa+ fa - m b )[( ai - b a 5 )(MCu 

S[{Ki + K3) m b\ 

+ faC\ 2 + m,gC )(aj + bj'j 5 )Q ij - G6Z]v(ki)(ki,k 3 ,m g ,m b .,m b .) , (A10) 

5M i 9 = ~ui , S , 2l "( fc 2)[K ~ 6,7 5 )(~ MC11- faC 12 + m- g C )( aj + brf)Qii 

s[(k 2 + k 3 y - m%\ 

-GSZ](fa+ fa + mbh x v(k 1 )Fx^eV(p 1 )e"(p 2 ){-k 2 ,-k 3 ,m g ,m bt ,m~ b .) , (All) 

6M " = ^mr^Mr 1 2 M k 2 )i\fa+ A3 - m b )[{ ai - b a 5 )(fac 11 + foc 12 + m- g c ) 

S[(K 2 + fc 3 j m b\ 

x(aj ■ + b j ^)Q ij -G5Z]v(k 1 )Fx^{p 1 )e v {p2)(k 1 ,k3,m g ,m b .,m b .) , (A12) 

5M ° 9 - [( fc2+P2 ) 2 - W D][(pT- fcl )-^] ^ 2h - (/;2 - h + mb)[{a > - b ^ ){MC " 

- hCu+ faC 12 + m g Co)(aj + bjj^Qij - GSZ] 

x(i»i- fa + m b Yi l ,v(k 1 )£ ll {p 1 )£ v (p 2 ){ki - p 1 ,k3,m g ,m b .,m b .) , (A13) 
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+m g C ){a j + b j j 5 )Q lj - GSZ](fa+ fa + m b )^ v 

x(i»i- fa + m b )7^v(k 1 )e^(p 1 )e u (p 2 )(-k2,-k 3 ,m g ,m ii ,m i .) , (A14) 

P2 + Pi) 2 - m 2 b )][{ P2 - fci) - m 2 ] 
- foC n + faC 12 + m g C )(a j + 6j7 5 )Qij - G5Z](p 2 - fa 

+m b )^ v v(k 1 )e ll (p 1 )e v (p 2 )(k 1 - p 2 ,k 3 ,m g ,m b ,m b .) , (A15) 



6M§» = ^ , 2 ^MrZ u,_ m *^ Yl^h- h + m b )[( ai - b^ifaCn 



6M!° = p 2 + fc3)2 ^g:_ fcl) _ m ^ (^)[K - ^7 5 )(- fcCn- AA, + m s C ) 
x(dj + bfi b )Q i0 - GSZ](fa+ fa + TOhh^J^- fa 

+m b )j„v(k 1 )e^(p 1 )s u (p 2 )(-k 2 ,-k 3 ,m g ,m bi ,m b .) , (A16) 

sm% 9 = — — — 2-wnrz 1 \ ^(^h^M- Pi+m b )-) v {fa+ fa 

[{k 2 - Pi) 2 - ml)\[{ki + k 3 ) - ml\ 

-m b )[(a,i - b l j 5 )(faC 11 + faC 12 + ■m g C )(a j + brf)Q i:j 

-G6Z}v(k 1 )e^(p 1 )e v (p2)(k 1 ,k s ,m- g ,m- b .,m~ b .) , (A17) 
— iGsT b T a 

SM * 9 = [(k 2 P 2 f - mWfaJ k 3 ) m 2f (^hAfa- h + m b h,(fa+ fa 
-m b )[(ai - b l j 5 )(faC 11 + faC 12 + ■m g C )(a j + &j7 5 )<2u 

-G8Z]v{k 1 )e^{p 1 )e v {p 2 ){k u k 3 ,m- g ,m b .,m bi ) , (A18) 

where we defined 

01,2 - (sin 6 b T cos b )/V2 , (A19) 
6 1)2 = (cos 6 b ± sin 6 b )/V2 , (A20) 
F\nv = (pi - P2)\9nv + (P2 + h + k 2 + k 3 )^x - (pi + ki + k 2 + k 3 ) v g\ ll . (A21) 

In the above, p\ andj?2 are respectively the momenta of the incoming quark and antiquark, ki,k 2 and £3 are respectively 
the outgoing b quark, b quark, and the Higgs boson momenta, s is the Mandelstam variable defined by s — (fa + 
k 2 + k 3 ) 2 = (pi + p 2 ) 2 , T a are the SU(3) color matrices, f abc are anti-symmetric structure constants of SU(3) color 
matrices, Gs = gg\Cp /16ir 2 with Cf — 4/3, and C and dj are the 3-point Feynman integrals [15] with their 
functional dependence indicated in the brackets following them. 
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